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ABSTRACT 

We investigate the excitation mechanisms of near-infrared [Fe n] and H2 emission lines 
observed in Active Galactic Nuclei (AGNs). We built a photoionization model grid consider- 
ing a two-component continuum, one accounts for the Big Bump component peaking at lRyd 
and another represents the X-ray source that dominates the continuum emission at high en- 
ergies. Photoionization models considering as ionizing source a spectral energy distribution 
obtained from photometric data of the Sy 2 Mrk 1066 taken from the literature were consid- 
ered. Results of these models were compared with a large sample of observational long-slit 
and Integral field Unit (IFU) spectroscopy data of the nuclear region for a sample of active 
objects. We found that the correlation between the observational [Feii]/ll.2570yum/Pay6 vs. 
H2/12.1218 yum/Bry is well reproduced by our models as well as the relationships that involve 
the H2 emission line ratios observed in the spectroscopic data. We conclude that the heating by 
X-rays produced by active nuclei can be considered a common and very important mechanism 
of excitation of [Fe n] and H2. 
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1 INTRODUCTION 

The excitation of the Narrow Line Region of Seyfert (Sy) galax- 
ies can reveal how radiation and mass outflows from the nu- 
cleus interact with circumnuclear gas. In particular, near-infrared 
(hereafter near-IR) observations are a powerful tool to investi- 
gate this issue, because the obscuration -which can affect the op- 
tical morphology of t he emitting gas region- is less importan t 
at these wavelengths dMulchaev et al .11 19961 : iFerruit et ail hoOO). 
Relevant emission-lines in the near-IR include [Fen] A 1.2570/jm 
and /t 1.6440 /mi, Hi lines such as Pa/3, and H 2 at A 1.9576 yum, 
A2.1218 fim, and /t2.3085/um, whic h can be used to map the 
gas ki n ematics and exci tation (e.g. Riffel & Storchi-Bergmann 
l2011bl : iRiffel et alj [2010). Nevertheless, the dominant excita- 
tion mechanisms of the [Fen] and H2 emission lines in the 
central regions of active galaxies are still unclear and have 
been the subject o f several recent studie s (e.g . Riffel et al 
l2010l 120081, 120061 : IStorchi-Bergmann et al.1 |2009l : Ifficks et al 



20091: iMtiller Sanchez etai 



Rodriguez-Ardila et al. 2005 



20091 : iRamos Almeida etai] 120091: 
l2004l : lDavies et alj|2007l) . ~ 



The H2 can be excited by two mechanisms: (i) fluorescent 
excitation through absorption of soft-UV photons (912-1 108 A) 
in the Lyman and Werner bands, existing both in star-forming 
regions and surrounding th e Active Galactic Nuclei (AGNs) 
(Black & van Dishoeck 1987) and (ii) collisional excitation due 
to the heating of the gas by sh ocks, the interaction of a r adio 
jet with the interstellar medium (Hollenbach & McKee 19891). o r 



by X-ray photons from the central AGN (Ma lonev et al. 11996). 
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Several studies based on intensity-line ratios teiffel et alJ) 20T 
IStorchi-Bergmann et al.ll2009l : lRodrfguez-Ardila et al.ll2005l,l200- 
have shown that collisional excitation processes dominate the H2 
emission surrounding AGN s. However, which is the dominant pro- 
cess is an open question. IVeilleux etai] ( 1 19971) . using J and K-band 
spectra of a sample of 33 Sy 2 galaxies, found that shocks associ- 
ated with nuclear outflows are a likely source of both [Fe 11] and 
H2 emission r ather than circumnuclear starbursts, as suggested by 
IOuillenetal] ( ll999l) . 

For the [Fen] emission, the [Fen] A 1.2570/jm/Pa/3 line ratio 
is generally used to investigate the main mechanism of excitation. 
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The value of this line ratio is controlled by the quotient of the vol- 
umes of partially and fully ionized gas regions , with [Fen] emis- 
sion being excited in the partially i onized gas (Mouri et al. 1990, 
1993l;lRodn'guez-Ardilaetai1 | 2005l; iRiffel et aljEoidl2008l [20061 ; 



Storchi-Bergmann et al. 20091). Such zones i n AGNs are created 



by X-ray emission (e.g. ISimpsonetaii ri996) and/or shock heat- 
ing of the gas by mass outflo ws from the nuclei w hich interact 
with the ambient clo uds (e.g. [F orbes & Ward 1993). This prob- 
lem was addressed by Mouri et al. (2000), who compared the val- 
ues of the line ratios [Fe n] A 1.2570 fim/Paj3 and [O i] A 6300 A/HJ3 
predicted by models, considering photoionization and shock heat- 
ing, with those observed in a sample of AGNs and Starburst 
galaxies. These authors pointed out that in AGNs, X-ray heat- 
ing is the most important [Fe n] excitation mechanism. However, 
Rodriguez-Ardil a"et al.l d2004l) . using near-IR spectroscopy of a 
sample of galaxies obtained with the Infrared Telescope Facility, 
found that X-ray excitation is enough to explain the H2 emis- 
sion and part of the [Fe 11] emission observed in Sy 1 galaxies, 
but fails to explain the emission of these elements in Sy2. For 
these objects, a combination of shocks and circumnuclear star- 
formation is required to explain these emissions. Moreover, it is 
not cle ar whether the [Fe n] an d H2 a re excited by the same mech- 
anism. Rodrigu ez-Ardila et al.l d2004l) found a correlation between 
the [Fe 11] A 1 .2570 //m/Pa/3 and the H 2 A 2. 1218 yum/Bry ratios, indi- 
cating that both sets of lines may be originated by a single dominant 
mechanism. However, high spatial resolution spectroscopy data 
from Integral Field Unit (IFU) of active galaxies indicate that the 
H2 and the [Fe 11] emitting gas have distinct flux distributions and 
kinematics, with the former being considered a tracer of the feed- 
ing of the AGN and the latter a tracer of its feedback dRiffel et al.l 



2010ll2009ll2008l;IStorchi-Bergmann et alj2009l ; lH"icks et al^Oo!* 
Mtille r Sanchez et alj |2009). This result indicates that the lines of 



these elements are formed in distinct regions. Although several 
works have investigated the excitation origins of H2 and the [Fe 11], 
it is still unknown whether a common mechanism can excite these 
elements. Fortunately, a large number of near-IR data of AGNs 
are currently available in the literature, which enables an extensive 
comparison with models yielding a more reliable conclusion about 
the likely dominant excitation mechanism of these emission lines. 

In this paper, we combined near-IR data of Sy galaxies ob- 
tained with IFU and long-slit spectroscopy with photoionization 
models to investigate the origin of the H 2 and [Fen]. In Section [2] 
we describe the observational data used in the analysis. The mod- 
elling procedures are presented in Sect.[3] In Sect. [4] the diagnostic 
diagrams used to compare the observational data with our model 
predictions are described. Results and discussion are presented in 
Sects. [5] and [6] respectively. A conclusion of the outcome is given 
in Sect. |7] 



2 OBSERVATIONAL DATA 

We compiled from the literature observational data of the nuclear 
region of active galaxies in the near-IR and optical spectral range 
obtained with long-slit and IFU spectroscopy. The selection cri- 
terion was the presence of bright infrared emission lines in their 
spectra. These data are described below. 



2.1 Long-slit data 

Near-IR e mission line intensity ratios of 35 active galaxies were ob - 
tained from Rodriguez -Ardila et al. l d2004|) . lReunanen et alJd2002h . 



lKnopetalJd200~u) . and lRiffel et al. (2006). This sample comprises 
long-slit data of 13 Sy 1 and 21 Sy 2 galaxies, along with 1 Quasar. 
The intensities of the near-IR [Fe 11] and H 2 emission lines observed 
in these objects were compared with our photoionization models. 
We also used the [O m] A 5007 A/H/3 and [0 1] A 6300 A/Ha line in- 
tensity ratios of about 600 000 emission-line galaxies listed in the 
MPA/JHU Data catalogue of the Sloan Digital Sky Survey DR7 re- 
lease (available at http://www.mpa-garching.mpg.de/SDSS/DR7/ 1. 

2.2 IFU data 

For this study we selected two Sy 1 galaxies, Mrkll57 and 
NGC4151, and two Sy2 galaxies (ES0 428-G14 and Mrk 1066). 
All of them were previously observed by our group using the IFU 
spectrographs of the Gemini telescopes. We selected these ob- 
jects because both J- and K-band spectroscopic data are available. 
The observations of Mrk 1066, Mrk 1157, and NGC4151 were 
performed using the Near-IR Integral field Spectrograph (NIFS; 
iMcGregor et alj|2003h on Gemini North, while ES0 428-G14 was 
observed with th e Gemini Near Infra-Red Spectrograph (GNIRS; 
lEliasetal .11 19981) on Gemini South. 



3 PHOTOIONIZATION MODEL 

To analyse the [Fen] and H2 excitation mechanisms, we built a grid 
of mo dels using the photoionization code Cloudy/08 dFerland et al.l 
1998), and then we compared the line intensity ratios predicted by 
them with those observed. The spectral energy distribution (SED) 
of the ionizing source used as input for the Cloudy code was a two- 
component continuum ranging from ~ 10 15 Hz to ~ 10 21 Hz. The 
shape of this SED is similar to the one observed in typical AGNs 
for that range. The first is the Big Bump component peaking at 
1 Ryd with a high-energy and an infrared exponential cutoff and 
the second one represents the X-ray source that dominates at high 
energies and is characterized by a power law with an index a x = 
-1. Its normalization was computed to produce the required value 
of the optical to X-ray spectral index a o x . This index describ es the 
continuum between 2 keV and 2500 A dZamoranietaI]|l98ll) . We 
assumed the default value of the Cloudy code a ox = - 1 .4, because 
that is about the average of the observed values, which are between 
-1.0 and -2.0, for the entire range of obser ved luminosities of AGNs 
dMiller et alj201ll;IZamorani et aljl98lh . 

The cosmic ray emission was considered in the models as a 
second ionizing source. Cosmic rays heat the ionized gas and pro- 
duce secondary ionizations in the neutral gas, which mostly in- 
crease the intensities of the H2 emission lines. We assumed a value 
of the H 2 ionization ra t e of 1 0~ 15 s" 1 , which is about the same rate 
found bv lMcCallet"ai]d2003h for a galactic line of sight. It is worth 
noting that the value of the cosmic ray rate must be estimated ob- 
ject by object. For example, ISuchkov et al.l d 19931) found for M 82 
a cosmic ray rate several times larger than the one in th e Milk Way. 
Gamm a ray observations of the starburst NGC 253 bv lAcero et al.l 
(2009) indicate a cosmic ray rate three orders of magnitude larger 
than that for the Milky Way. Also, molecular data of star forming 
galaxies, such as Arp 220, show evidence for extremely high cos- 
mic ray rates yielded b y the UV emission from supernova remnants 
(Meijerink e t alj201ll) . 

We computed a sequence of models assuming an electron 
density N s = 10 4 cm 4 , ionization parameter U in the range 
-4.0 < logf/ < -1.0 defined as U = Q io „l4nR\nc, where Q ion 
is the number of hydrogen ionizing photons emitted per second by 
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Table 1. Fe/O and O/H gas phase abundances assumed in the models. 



Metallicity (Z/Z G ) 12+log(0/H) log(Fe/0) 



2 9.0 -1.47(al) -1.94(a2) -2.24(a3) 

1 8.69 -1.77(bl) -2.24(b2) -2.77(b3) 

0.5 8.38 -2.15(d) -2.54(c2) -2.76(c3) 



the ionizing source, R$ is the Stromgren radius (in cm), n is the 
particle density (in crrr 3 ), and c is the speed of light. The cho- 
sen range of thes e values for U is typical of narrow-line regions of 
Sy galaxies (e.g. iFerland & Netzeiill983h . The Ht emission lines 
are very dependent on the electron density value assumed in the 
models. For example, when N e varies from 10 3 to 10 5 cm , the 
logarithm of the H 2 A 2.1218yum/Bry emission line intensity ratio 
span about 2.6 dex. The value N e = 10 4 cnT 3 , assumed in our mod - 
els, is a mean value from those considered bv lMourfe t al. (200G). 
We considered in our models three values of 12+log(0/H)= 8.38, 
8.69, and 9.00, which correspond to values o f the metallicity 0.5, 
1, and 2 times the solar value published by lAllende Prieto et all 
d200lh . The abundances of other metals in the nebula were scaled 
linearly to the solar metal composition through the comparison 
of the oxygen abundances, with the exception of the N and Fe 
abundances. The nitrogen abundance was taken from the relation 
log(N/O)=log(0.034+120O/H) of Vila- Costas & Edmunds (1993). 
The Fe/O abunda nce ratio has a large scatter for a fixed O/H value 
dlzotov et al.ll2006h and its value is uncertain because the Fe and 
O ab undances in grains are poorly known dPeimbert & Peimbertl 
l20ld) . Thus, we varied the Fe/O abundance ratio by about 0.7 dex 
on each metallicity. 

The presence of internal dust was considered and the grain 
abundances (van Hoof et al. 2001) were also linearly scaled with 
the oxygen abundance. To take into account the depletion of re- 
fractory elements onto dust grains, the abundances of Mg, Al, Ca, 
Ni, and Na were reduced by a factor of 10, and Si by a factor of 2 
relative to adopte d abundances in each model in accordance with 
iGarnett et ail Jl995h . In Table[Q the O/H and Fe/O abundance val- 
ues of the gas phase assumed in the models are sh own. The model 
of the H 2 molecule described by Shaw et al. ( 2005) and the model 
of the Fe + ion described by IVerner et alj ([1999), which consider 
371 energy levels, were assumed in our computations. The outer 
radius of the modelled nebula is that where the temperature falls 
below 1000 K. 
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Figure 1. Diagnostic diagram showing the observational data taken from 
the literature (see Sect.|2j and results from the grid of photoionization mod- 
els (see Sect|3j. Solid lines connect curves of iso-Z, while dotted lines con- 
nect curves of iso-l/. The values of log U and Z are indicated. The three 
different lines for each Z correspond to the different assumed values of the 
Fe/O as indicated by the labels (see Table[TJ. Circles, squares, and star rep- 
resent Syl, Sy2, and quasar data, respectively. The typical error bar (not 
shown) of the emission line ratios is about 10 %. 

H 2 A2.247 y um/i2.121 y um (Fig. Hi — llvlouril d 19941) proposed 
these diagrams to separate gas emission yielded by shocks from 
emission caused by fluorescence. The drawback in using the 
H 2 /l 1.957yum//l 2.121 yum ratio is that the H 2 /1 1.957yum may be 
affected by telluric bands of H 2 Q and CQ 2 , or blended wi th the 
[Sirv] A 1.963 fj.m emission line 1 Rodrfguez-Ardil a et al.ll2005l) . 

• [O m] A 5007A/H/3 vs. [O i] A 6300A/Htt (Fig. [3} — This di- 
agram was suggested bv lBaldwinetalJdl98ll) to separate objects 
according to their primary excitation mechanisms, i.e. (a) photoion- 
ization by stars, (b) photoionization by a power law continuum 
source or (c) shock heating. In particular, the [O i] A 6300 AfHa 
line ratio is greatly increased b y the presence of s hocking gas, even 
when it has low velocities (e.g. lAllen et al.l l2008). 



4 DIAGNOSTIC DIAGRAMS 

We used four diagnostic diagrams containing predicted and ob- 
served emission line ratios of the [Fe n] , H 2 , [O in] , and [O i] which 
are described below. 

• [Fen] Al. 2570 //m/Pa/? vs. H 2 /l 2.1218^m/B r- v (Fi g. [Q 
— Diagnostic diagram sugge sted by lLarkin et al.l J 19981) and 
Rodriguez- Ardil a"et al.1 d2004l) to separate galaxi es accord- 
ing to their level of nuclear activity. Recently, iRifFel et all 
j2010l) constructed this diagram with spatially resolved 
IFU data of an AGN. Typical values for the nucleus of 
Sy galaxies are 0.6 < [Fen] A 1.2570 nm/Pafi< 2 an d 
0.6 <H 2 ^2.1218^m/Bry£ 2.0 jRodriguez-Ardila et alll2005h . 
This [Fe n]/Pa/? is very dependent on the Fe/O abundance while the 
H 2 emission lines are dependent on the ionization parameter. 

• H 2 ^ 1.957 /zm//l2.121 y um and U 2 A 2.033 fim/A 2.223 /mi vs. 



5 RESULTS 

5.1 Integrated spectra 

In Figs. Q] and [2] we show the first three diagnostic diagrams de- 
scribed above containing the results of our grid of photoionization 
models and the data sample. Sy 1, Sy 2 and quasar are represented 
by different symbols. For the IFU data, the emission line ratios rep- 
resented in these Figs, were estimated by integrating the spaxels 
inside a central aperture of 0.5"x0.5" for each galaxy, with excep- 
tion of ES0428-G14 for which an aperture of 0.75"x0.75" was 
considered. These values are presented in Table|2] 

In Fig. Q] we can see that almost all the observational ratios 
are within the parameter space defined by our grid of photoioniza- 
tion models. A lower metallicity than those assumed in our mod- 
els is required to reproduce the data of the galaxies out of the 
grid. The observed correlation between [Fen] A 1.2570yum/Pa/? and 
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H 2 2.247/im/2.121/im 

Figure 2. As in Fig. [T] for H2 emission lines. The arrow indi- 
cates the direction in which the ionization parameter increases. Circles, 
squares, and star represent Syl, Sy2, and quasar data, respectively. The 
hatched area represen t s the region occupied by shock model results from 
iHollenbach & McKed fl98Sj) . 



Table 2. Integrated line ratio intensities of IFU data 



Object 


[FenU1.2570pm/Pa/3 


H 2 /12.1218^m/Bry 


ES0428-G14 


0.75 


1.10 


Mrk 1066 


0.52 


0.96 


Mrkl 157 


0.73 


2.24 


NGC4151 


0.45 


0.26 



H2 A 2.1218/im/Bry is explained by an increase in metallicity and 
ionization parameter. Noteworthy that the parameter space defined 
by the models built using Z = 0.5 Z Q is almost completely con- 
tained in the one denned by the models built using the solar metal- 
licity. 

In the case of the diagnostic diagrams that only involve H2 
line ratios (Fig. [2}, the photoionization models are slightly de- 
pendent on the assumed metallicities, covering almost the same 
parameter space, and strongly dependent on variations in the 
ionization parameter. Taking into account the observational er- 
ror bars, our models are in good agreement with the observed 
H 2 /l 2.033 /jm//i 2.223 fim and H 2 A 2.247 fjm/A 2.121 /jm ratios (up- 
per panel). On the other hand, in the lower panel of Fig. [2] can be 
noticed larger dispersion of the observational data which is not well 
reproduce by the models. This dispersion could be the result of a 
contamination of the measurements of the H2 A 1.957 jim emission 
line intensities due to a blend with the [Siiv] A 1.963 yum line (as 
explained above). Therefore, the predicted H 2 /11.957yum intensi- 
ties are somewhat lower than the observed ones. In Fig. [2] we also 
show the area occupied by the theoretical intensities of the line 
ratio H 2 /i2.247^m//12.121 from shock models performed by 
IHollenbach & McKee Jl989h . These authors computed emission- 
line spectra of steady interstellar shocks in molecular gas con- 
sidering velocities from 30 to 150 km/s and particle densities of 
10 3 - 10 6 ctrr 3 . We can see that most of the objects of our sample 
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Figure 3. [O m A 5007 A/H/3 vs. [0 1 A 6300 A/Ha diagnostic diagram. The 
yellow line sep arates objects ion ized by massive stars from those containing 
active nucleus I Kewley et al. 2001). Blue, green and red solid lines are as in 
FigU Points represent emission-line galaxies listed in the MPA/JHU Data 
catalogue of the Sloan Digital Sky Survey DR7 release (see Sect.[2). 

appear to have the X-rays as main ionizing source while for the re- 
maining ones a composite ionization by X-rays and shock can be 
considered. 

Fig.[3]shows the [O in] A 5007 A/H/3 vs. [0 1 A 6300 A/Ha diag- 
nostic diagram. In this Fig. we can see that the observational data of 
AGNs are well describe by our models. If our models use the lower 
values of the ionization parameter (log U < -3.5; these models 
are not shown), we can extend the parameter space to include the 
objects that have values of the logarithm of the [O in] A 5007 A/H/3 
ratio lower than zero. As in the case of the [Fen] A 1.2570yum/Pa/3 
and H 2 A 2.1218yum/Bry diagnostic diagram (Fig.QJ, the parameter 
space of the models with solar metallicity almost contains that of 
the Z = 0.5 Z Q models. 

5.2 IFU data 

We plot the [Fe 11] A 1 .2570 //m/Pa/3 and H 2 A 2. 1 2 1 8 yum/Bry diag- 
nostic diagram for each spaxel of our four objects with our model 
results (see upper panels of Figs. [4] and [5j. In these Figs., the 
spaxel data are separated by their ionization mechanism accord- 
ing to the place in the diagnostic diagram, with different colours 
for each mechanism. The different ionization mechanism zones 
are delimited in the Figs, b y dashed-lines, following the work of 
iRodnguez-Ardila et al.l fl2004h . The spaxels showing typical values 
of starbursts, Seyferts, and low-ionization nuclear emission-line re- 
gions (LINERs) are represented by green open circles, black filled 
circles, and red crosses, respectively. With the same colour code, 
we show the spatial position of each spaxel in the IFU field of view 
(see lower panels of Figs. [4] and [5). Our models completely repre- 
sent the region occupied by Seyfert and LINERs data. 

6 DISCUSSION 

The excitation mechanism of the near-IR emission lines of the 
[Fen] and H 2 in active galaxies have been the subject of sev- 
eral works. For example, iMouri et"al] d2000l) compared results of 
models considering photoionization by X-rays and shock heating 
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Figure 4. Top panels: [Fe 11] A 1.2570 /im/Pa/3 vs. H 2 A 2. 1212 ji/m/Bry line-ratio diagnostic diagram for Mrk 1157 (left) and ES0428-G14 (right). The dashed 
lines delimit regions with ratios typical of Starbursts (green open circles), Seyferts (black filled circles) and LINERs (red crosses). Blue, green and red solid 
lines are as in Fig.[T] Bottom panels: spatial position of each spaxel in the IFU field of view from the diagnostic diagram. 



with observational data of AGNs and starburst galaxies. These 
authors built their models considering large ranges in shock ve- 
locities, gas density, metallicity, and different ionizing continua. 
Mouri and collaborators showed that the [Fen] emission is en- 
hanced when a partially ionized zone is produced by photoion- 
ization by X-rays (described by a power-law) and shock heating. 
These two processes can be discriminated by the electron tem- 
perature of the [Fen] region: 8000 K in heating by X-rays and 
6000 K in shock heating. Com paring the e l ectron temperature of 
the [Fen] region estimated by iThompsonl d 19951) for NGC4151 
(8000 < T c < 12000 K) with their models. iMouri et all fcOOOh 
showed that, at least for this galaxy, it indicates that X-rays are 
the more important mechanism to yield the [Fe n] flux. These au- 
thors arrived to the same conclusion using the [O i] A 6300 A/Ha vs. 

/Pa/3 diagnostic diagram. A similar result was 
Jackson & BeswickI (2007) by analysfng J-band 



[Fen]/11.2570yUM 

ALSO OBTAINED BY 
SPECTRA OF THREE Sy 2 GALAXIES. 



For our models, we assumed an incident continuum whose 
shape is given by two components, a Big Bump and an X-ray power 
law, varying the Fe/O abundance. With these models that do not 
consider shock heating, we are able to explain the observational 
data. Nevertheless, we do not exclude some contribution by shock 
heating to the [Fen] emission. Comparing our models with the 
SDSS DR7 emission-line galaxies (Fig. [3}, we are able to describe 
the [O i][Ha line ratio observed in AGNs. This diagram cannot be 
used to discriminate [Fe n] excitation mechanism, nevertheless, we 
must take into account that the [O i]/Ha line ratio is shock sensitive. 
Hence, although shock contribution in the ionization of Fe cannot 
be ruled out, models considering a continuum described by a Big 
Bump and an X-ray power law as the ionization source can also re- 
produce the [Fen] emission lines as well as the behaviour of shock 
sensitive emission lines such as [O i] A 6300 A. Analysing I FU ob- 
servations of the Sy galaxy NGC 4151. iTurner et all d2002l) found 
that the [Fe n] emission mainly arises in the visible narrow-line re- 
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Figure 5. As in Fig.g]but for Mrk 1066 (left) and NGC4151 (right). 



gion in which the dominant excitation mechanism is th e photoion- 
ization by collimated X-ray emission from the nucleus. lOliva et al.l 
feooih pointed out that in regions where shocks are the dominant 
mechanism the iron-based grains are destroyed but the phosphorus 
is not, yielding a larger [Fen] A 1.2570yum/[Pn] A 1.188yum line- 
ratio intensity than that observed in the regions dominated by X- 
ray. In order to verify this, in Fig.|6]we show a histogram containing 
this observed line i ntensity ratio for 17 Seyfert galaxies, 5 Sy 1 an d 
12 Sy 2, taken fromljackson & Beswickl ( l2007h jRiffel eTallhoOfj) . 
and lOliva et a I] j200lh . It can be seen that the [Fen]/[Pn] for Sy 
galaxies ranges from 1.5 to 6 (with a mean value of 2.7). The mean 
value of this ratio is about 20 for SNRs, which indic ate that the 
emitt ing gas has recently passed through a fast shock jOliva et all 
1200 lb . Therefore, these results confirm that shocks have little influ- 
ence on the [Fe n] emission. 

Regarding the H2, this molecule can be excited via three dis- 
tinct mechanisms: (1) UV fluorescence, where photons with A > 
912 A are absorbed by the H 2 molecule and then re-emitted, result- 
ing in the population of various vibro-rotational levels, (2) shocks, 



where high-velocity gas motions heat and accelerate this molecule; 
and (3) X-ray illumination, where hard X-ray photons penetrate 
deep into molecular clouds, heating large amounts of molecular 
gas resulting in the IT emission (see lRodrfguez-Ardila et al 1120041 
and references therein). Rodriguez- Ardila and collaborators used 
the diagrams shown in Fig.[2]to compare observational data of 22 
objects with models considering a thermal emission, a non-thermal 
UV excitation, a thermal UV excitation, and a mixture of thermal 
and low-density fluorescence. These authors found that for 4 ob- 
jects the excitation mechanism is clearly thermal, while for the re- 
maining objects a mixing with a non-thermal process cannot be 
discarded, even though the results point out to a dominant thermal 
mechanism. 

To analyse the relative weight of the X-ray emission with re- 
spect to the other model components (mainly with fluorescence and 
UV photons), not only for the IT emission but also for the [Fen], 
we made models fixing all parameters with the exception of the a ox 
value (see Fig. |7J, which is related to the X-ray power law nor- 
malization (see Section $3}. We assumed Z = Z Q and log U = -2.5 
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[Fell] 1 .257u.m/[PII] 1 .88^m 

Figure 6. Histogram showing the [Fe n]/[P n] emission line intensity ratios 
of a sample of objects collected from the literature. 



since the models built using the solar metallicity and this value of 
the ionization parameter cover almost all the parameter space occu- 
pied by the ob servational data (see Fig ,[T), Taking into account the 
a ox definition (Tanan baum et alJI 19791) . which fixes the Big Bump 
parameters, a decrement of the a ox value implies that the amount 
of the X-rays emitted by the source decreases. In Fig.UJwe can see 
that our models with a m = - 1 .4 reproduce well the observational 
AGN data. Nevertheless, when we use lower values of this parame- 
ter, the ratios predicted b y the models go out of the reg ion typically 
occupied by the AGNs JRodrfguez-Ardila et alj|2004l). Therefore , 
our models favour the scenario suggested bv lMalonev et al. I J 19961) . 
where the H 2 molecule emission is mainly governed by photons 
emitted at X-ray wavelengths from the central AGN. This also can 
be inferred from the dependence of the H2 emission lines on the 
ionization parameter U . 

To verify if shock models can fit the observational data, we 
compared shock model results by iHollenbach & McKed jl989h 
with our sample (see Fig. [2). Only few observational points are 
located in the area occupied by these shock models and, even in 
these cases, models considering X-rays also describe the data. 

On the other hand, varying in our models the H2 ioniza- 
tion rate by cosmic rays by a factor of 200, we found that the 
H 2 A 2.1218/jm/Bry line ratio only increases by about 0.15 dex, 
which shows that the additional ionization by cosmic rays has little 
influence on the H 2 emission lines. 

A simple scenario where both [Fen] and H 2 emissions are 
mainly due to the X-ray continuum coming from the ac tive nucleus 
has al so be en proposed by othe r authors. For example, iBlietz et"al] 
dl994h and iKnop et al.l ( feOOlh showed that X-rays from the nu- 
cleus can heat the gas located in the narrow line region driving 
the [Fe 11] and H 2 emission. Because 98 % of the iron is tied up 
in dust grains, this process must free the iro n through dust destruc- 
tion a nd yet not destroy the H 2 molecules (Rodriguez- Ardi la et al.l 
1 20041) . These authors computed the emergent [Fe n] A 1.2570 fim 
and H 2 A 2. 1 2 1 8 pm flux using the X-ray models bv lMalonev et al.l 
dl996h and compared their predictions with observational data of 
seven objects. They found that X-ray heating can only explain 
a fraction of the [Fen] and H 2 emission, and they stated that 
the discrepancy found can be alleviated if the emitting gas is lo- 
cated closer than the distance adopted in their models. The X- 
ray data, provided by the XMM-Newton an d Chandra space tele- 
scopes, and their detailed analysis (see e.g. IPiconcelli etai]|2005l ; 




ot(ox) 

Figure 7. Model results using solar metallicity, N e = 10 4 cirT 3 , logU = - 
2.5 and varying only the a ox parameter to see the influence of the X-rays on 
the [Fe 11] A 1.2570pm/Pa/3 and H 2 A 2.1218/jm/Bry emission line ratios. To 
delimit the region occupied by the AGNs we follow Rodn'guez-Ardila et al. 
120041) . 
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Figure 8. Spectral energy distribution at the Schwarzschild radio 
(10~ 5 pc) of the Sy 1 galaxy Mrk 1066 used as the photoionization 
source for some models of this g alaxy. We assumed a galaxy dis- 
tance of_50Mpc i Mould et al. 2000). The photometric data were taken 



from Dressel & Condon 1 1978), Moshir et al 



1 1990) , de Vaucouleurs et al 
. 119961). ICondon et all feOO: 
I Skrutskie et alj 12003|), iBraatz et all 120041), iGuainazzi et alj 
iMufioz-Marin et alj 120071) . and lcardamone et alj 120071) ' 



1 1991) 



Beck er et all Jl99lh- iDouglas et all 



Longinotti et al. l2007l : Bianchietal. 20091 : iKrongold et alj|2009t 
Carda cTet alj l201ll : ICorral et al.l 2011 , and references therein), 



provide information about the continuum shape and the particu- 
lar spectral features of the AGNs in this wavelength range. For 
Mrk 1066, we compared the results obtained using this simple sce- 
nario that only involves a continuum modelled by a Big Bump 
and an X-ray power law with those obtained using its intrinsic 
SED. We built the observational SED taking the photometric data 
from the NASA/IPAC Extragalactic Database (NED), following 
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log(H 2 A2.121/im/Br7) 

Figure 9. Comparison between the grid model results shown in Fig.[T](solid 
lines) and the models built considering the semi-empirical SED of Mrk 1066 
shown in Fig.[8](dashed lines). 



by iRiffel et alj ( bOOol I2008L |2010|) and IStorchi-Bergmann et al.l 
(2009), respectively. Howev er, these authors did not reach con- 
clusive results. For example, Riffel et al. ( 2006) suggested that the 
[Fen] excitation in ES0 428-G14 is mainly due to shocks. Never- 
theless, the detailed analysis performed in the present work con- 
fronting our models with the IFU data shows that X-rays are a 
more reliable dominant excitation mechanism even in the case of 
ESQ428-G14. 



7 CONCLUSIONS 

In this work we show that a photoionization model grid built by 
adopting a continuum source characterized by two components, 
one accounting for the Big Bump component peaking at 1 Ryd and 
the other describing the X-rays emission, is able to reproduce the 
[Fe n] and H 2 infrared emission lines of a sample of AGNs. Testing 
the influence of the X-rays on the intensity of these emission lines, 
we found that a decrement in the X-ray content of the continuum 
source translates into a weakening of these lines, and the models 
are no longer compatible with the observations. This implies that 
the heating by the X-ray emission from the active nuclei can be 
considered as the most important mechanism of excitation for the 
IR emission lines of these elements. 



ICardaci et al] (§009). To enhance the number of points of the SED 
as needed by Cloudy, we performed a linear interpolation among 
the semi-empirical points (see Fig.[8jl. We built a new grid of pho- 
toionization models under the same assumptions of abundances, 
ionization parameters and density, but only for one value of the 
Fe/O ratio for each metallicity. In Fig.[9]the predictions of our mod- 
els using the SED of Mrk 1066 and the model results presented in 
Fig. [2 assuming the same Fe/O abundance as the Mrk 1066 mod- 
els are shown. The model results derived using the two different 
ionizing sources are mostly in agreement. 

The semi-empirical SED of Mrk 1066 includes not only the 
range covered by the Cloudy model but also the radio and IR wave- 
lengths. Hence, the agreement between solid and dashed lines in 
Fig. |9] only indicates that the assumed multicomponent model is 
a good representation of the AGN continuum when studying the 
[Fe n] and Ht emission. 

Recent resolved integral field spectroscopy of the central 
region of active galaxies shows that the ionized (in particu- 
lar the [Fen] emitting gas) and the molecular (traced by the 
H2 emission) gas have distinct flux distributions and kinemat- 
ics. The molecular component is more restricted to the plane of 
the galaxies and the ionized one extends to high latitudes above 
it, which is in most cases co - spatia l ly with the radio jet (e.g., 
Riffel et alj|2006l. l2008l I2009L 1201 d: Riffel & Storchi-BergmanrJ 



201 lal fbl: IStorchi- Bergmann et al. 2009i |2010|) . Usually the [Fen] 



has an enhancement in flux and velocity dispersion in regions sur- 
rounding the radio structure, suggesting that the radio jet plays an 
important role in the [Fen] emission. Our models are able to re- 
produce the [Fen] emission of active galaxies without consider- 
ing shock excitation by the radio jet. Thus, the enhancement in the 
[Fe 11] flux in the vicinity of radio structures can be interpreted as 
being due to an enhancement in the gas density, caused by the in- 
teraction of the radio jet with the emitting gas, and mainly excited 
by X-rays from the central engine. 

The main exciting mechanism of infrared emission lines of 
ES0 428-G14, Mrk 1157, Mrk 1066 and NGC4151 was discussed 
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